POLYNOMIAL FIT FOR JITTER SEPARATION 



BACKGROUND OF THE INVENTION 
The present invention relates to jitter analysis. 

Characterization of the transient behavior of e.g. high-speed digital circuits, i.e. the 
5 transition of a digital signal from one logical state / level to another logical state / level, 
has become growingly important for designing as well as manufacturing such digital 
circuits. With a continuously increasing level of SOC (system on a chip) integration 
multiple transistors switch at the same time thereby disadvantageously influencing 
timing characteristics of the signals generated by means of switching noise. Further, 
10 electrical fields generated by nearby clock distribution lines, etc. can also cause timing 
instabilities such as jitter. Accordingly, an amount of jitter that can be tolerated or that 
characterizes the transmission quality, e.g., of an lO-device is used as a figure of merit 
in specifications or standardizations. 

ITU-T G.701 defines jitter as short-term non-cumulative variations of the significant 
15 instants of a digital signal from their ideal positions in time. The significant instant can 
be any convenient, easily identifiable point on the signal such as the rising or falling 
edge of a pulse or the sampling instant. By plotting the relative displacement of the 
instants within a real pulse train that has some timing jitter with respect to an ideal 
pulse train, the so-called jitter time function can be obtained. Alternatively, the jitter 
20 spectrum can be displayed in the frequency domain in order to provide jitter analysis. 
In a probabilistic approach jitter-histograms showing the likelihood for a transition can 
be established and analyzed. 

In order to characterize devices with respect to jitter semiconductor manufacturers 
expect to use measurement equipment that is able to separate different types of jitter 
25 in order to allow an identification of sources of error and to provide further 
improvements of product quality. In particular, it is generally desirable to separate 
random jitter (RJ) from deterministic jitter (DJ). However, existing techniques such as 
that disclosed, e.g., in Li, M P., Wilstrup, J., Jessen, R., Petrich, D. in: "A New Method 
for Jitter Decomposition through its Distribution Tail Fitting", Proc. Intl. Test Conf., pp. 
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788 - 794 (1999) are very time consuming, need many measurements to be taken per 
data point and therefore do not meet the requirements of high volume production. 

SUMMARY OF THE INVENTION 

It is an object of the invention to provide improved jitter separation. The object is solved 
5 by the independent claims. Preferred embodiments are shown by the dependent 
claims. 

According to the present invention, a bit error rate function being obtained by means 
of measurement using, e.g., a commonly available bit error rate tester is analyzed to 
separate and determine the amount each of deterministic and random jitter. For this 
10 purpose a plurality of bit error rate values are first measured as a function of compare 
time in a bit time interval. Successive timing points are selected each therefore. The 
term "plurality" may include the number of two, three or a larger count of timing points. 

As explained in detail, e.g. in NCITS T11.2 / Project 1230 "Fibre Channel - 
Methodologies for Jitter and Signal Quality Specification" , Rev. 5.0, Feb. 12, 2002, a 

15 bit error rate (BER) represents a cumulative probability of a transition causing an error 
due to jitter. To determine a BER, a sample point in time is usually swept along the bit 
time axis with each measurement collecting the events of a transition of a digital signal 
level prior to the center of the data eye and the results are then compared with 
expected values. This is valid if one tests for the decreasing data eye from the "left 

20 side" of a bit window. The opposite side of the bit time interval can similarly be tested 
by collecting any transition event occurring after the current timing point under 
investigation until a sufficient time after the bit window has ended and then - again - 
comparing said results. 

Once the bit error rate values are obtained for the respective timing points a polynomial 
25 function is fit to the data obtained on each side of the data window. A perception 
underlies the invention, that due to its random jitter portion particularly the cumulative 
bit error rate function may be well represented by a polynomial function, e.g., of limited 
order not only in its deep tail, but also at moderate time distances from an expected, 
ideal transition. 
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Thus, according to one aspect of the invention, measuring bit error values at such 
timing points, which will provide BER values lying substantially within the range of 10" 4 
and 10" 6 , comparatively short measurement times become feasible in contrast to 
hitherto existing fitting methods, in which there were to measure BER values down to 
5 10~ 12 reflecting rare occurrences of transition events in order to, e.g., obtain sufficiently 
reliable results for the individual jitter histogram representations. 

Random jitter is assumed to cause a Gaussian type probability distribution of transition 
events that is unbounded in time and being centered at the expected transition, if 
deterministic jitter is absent. The amount of random jitter is related to the standard 

10 deviation of the distribution function. However, in case that deterministic jitter is 
present, the distribution function splits into two separated Gaussian peaks by means 
of convolution of two Dirac peaks representing two bounded levels of deterministic 
events with said Gaussian type probability distribution. The corresponding BER 
characteristic is therefore equivalent to the integral of this convoluted function. The 

15 amount of deterministic is related to the distance of the two peaks. 

What is derived in the polynomial fit to the BER according to the present invention are 
polynomial coefficients. These are evaluated to identify the amount each of random 
jitter and deterministic jitter, separately. 

In a preferred embodiment the bit error rate values measured are transformed prior to 
20 being fitted with a polynomial function. In this transformation advantage is taken of the 
form of the cumulative distribution function, i.e. BER characteristic, which can be 
approximated by a linear regression for BER values lower than a threshold of, e.g., 10" 
4 , if at first the, e.g. decadic logarithm is taken of bit error rate values and at second a 
square root is taken of the values such derived. The slope and, e.g., the crossover of 
25 the resulting trend- or regression line with the ordinate axis can be analyzed and used 
to determine and separate random (slope, or first order polynomial) and deterministic 
jitter (ordinate axis crossover, or zeroth order polynomial), respectively. 

The invention can be partly or entirely embodied or supported by one or more suitable 
software programs, which can be stored on or otherwise provided by any kind of data 
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carrier, and which might be executed in or by any suitable data processing unit. 
Software programs or routines are preferably applied to processing units of commonly 
available Bit Error Rate Testers (BERT). As a consequence of the short measurement 
times provided by the present invention in situ output of jitter analysis results is enabled 
rather than spending efforts in performing a costly post-processing analysis, separately. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Other objects and many of the attendant advantages of the present invention will be 
readily appreciated and become better understood by reference to the following 
detailed description when considering in connection with the accompanied drawings. 
Features that are substantially or functionally equal or similar will be referred to with the 
same reference signs. 

Fig. 1 shows a probability distribution histogram for transitions of a logical level of 
a digital signal with respect to an ideal timing event; 

Fig. 2 shows a cumulative BER function of the histogram curve as measured (thin 
line) being displayed in Fig. 1 in logarithmic representation, and a 
polynomial fit to the data (bold line); 

Fig. 3 shows the cumulative BER function as displayed in Fig. 2, but transformed 
by changing sign and square rooting the data as measured (thin line) and 
as those fitted (bold line); 

Fig. 4 shows the cumulative BER function as displayed in Fig. 2, but a polynomial 
fit (bold line) performed considering an offset of the baseline; 

Fig. 5 shows the cumulative BER function as displayed in Fig. 2, but compared with 
actually measured data of three timing points used as input for the 
polynomial fit (bold line) for determining total jitter; 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS OF THE INVENTION 

For simplicity and clarity reasons, in what follows theoretical data representations have 
been chosen to demonstrate the present invention. However, it is clear to the person 



skilled in the art, that the invention can advantageously be applied to actually 
measured data distributions as well. In the drawings, the theoretical - or equivalently: 
measured - data representations are marked with an "m", while polynomial fits are 
marked with an T, respectively. 

5 Fig. 1 shows a theoretical probability distribution histogram for transitions of a logical 
level of a digital signal with respect to an ideal timing event vs. time units Ul, wherein 
1.0 Ul is the length of a bit window. Due to the presence of random as well as 
deterministic jitter a bimodal distribution that is approximately symmetric about the ideal 
timing event (t=0) occurs. Integrating this probability distribution function (PDF) with 
10 respect to time t, the cumulated distribution function - or bit error rate (BER) function 
can be obtained, which is displayed in Fig. 2. 

Since the probability density function of such a bimodal jitter is represented by two 
Dirac peaks the overall jitter distribution is a sum of two Gaussian PDF's resulting from 
the convolution of the Gaussian shape with the Dirac pulses: 

, t -fell , , 
4V2^-cr 1 4V2;r.<7 2 

wherein o-| = 02 is the standard deviation representing the random jitter, x refers to the 
time variable, and u<| , U2 are tne deviations from an ideal timing event resulting from 
deterministic jitter. In this model representation, typical values for ai=a 2 =0.02 Ul 
(random jitter), and ^=^=0.1 Ul and n 2 -Hi=0.2 Ul (deterministic jitter have been 
20 applied). 

The bit error rate values displayed as the thin line in Fig. 2 have been transformed into 
a logarithmic form, log-|Q (BER). As can be inferred from the functional behavior of the 
analytical expression of the PDF, a second order polynomial can easily be fitted to the 
logarithmized form of the BER for at least larger values of the time coordinate (bold 
25 line). 

In order to linearize the fit, a further transform is applied by square rooting the bit error 
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rate values, after changing the sign of the logarithm ized values. The result is shown in 
Fig. 3. The polynomial fit (bold line) is derived by using bit error rate values selected 
from a time interval between 0.18 Ul and 0.40 Ul - or a 0.08 Ul to 0.30 Ul time distance 
after the actual occurrence of the second Dirac peak that is due to deterministic jitter. 
5 Since the interval is to be adapted to the actual amount of random jitter, i.e. the value 
of said standard deviation o, an iterative procedure pre-setting a first reasonable guess 
of a standard deviation and then refining the interval appropriately. In the present case 
the interval starts at 4 a. 

Due to deterministic jitter the quadratic portion of the Gauss-profile is shifted to the right 
10 by the amount of m, which can be seen in Fig. 4. The vertical shift of the quadratic 
section is small, which is marked with an "O" in Fig. 4. When we deploy the technique 
of taking the square root of the graph in Fig. 2 we neglect in this embodiment a small 
vertical offset of this quadratic portion. However, on the one hand side, it has been 
found that this will introduce only negligible errors to the individual jitter determinations, 
15 on the other hand side a more detailed analysis regarding this offset we be given below 
in another embodiment. 

Therefore, a trendline T fit by linear regression is enabled and an approximation of the 
BER section assumed to show quadratic behaviour, i.e. larger than 0.08 Ul in the 
present example, to the logarithmic representation of the BER is carried out. An 
20 extraction of the amount of random jitter can simply be performed by determining the 
slope of the trendline T and then relating the slope to the argument of a complementary 
error function (erfc), which gives an analytical expression of one of the integrated 
Gaussian profiles, i.e. an idealized BER in the absence of deterministic jitter: 



Or equivalently, an assumption is made, that the contribution of the first probability 
peak to the tail of the second peak is negligible. Under this assumption, a relationship 
between the slope of the trendline T determined according to the fit and the standard 
deviation o can be established. In this embodiment, a o-value of 0.0198 Ul is 
30 recovered, which is in close agreement with the value of 0.02 input to the theoretical 




25 
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PDF. 



With respect to the determination of the amount of deterministic jitter the zeroth order 
coefficient of the trendline T is considered: in case of jitter only being represented by 
random jitter the line would intersect the ordinate axis at approximately the baseline B 
5 of the BER characteristic, i.e. at the value of sqrt(-log-|o(0.5)) = 0.54866 (Note that the 
probability of an event occurring at all is just 0.5 and not 1). If deterministic jitter DJ is 
present, the quadratic section of the BER is shifted, and the crossover points of the 
trendline T with the ordinate axis as well as the baseline B change. The latter crossover 
point gives a measure for deterministic jitter DJ on the time axis. This crossover point 
10 can easily be derived from the trendline T and then be identified with \x v Without 
restriction the opposite boundary of a data eye can be treated in the same way thus 
obtaining \x 2 . The result is graphically displayed in Fig. 3. In spite of the small offset 
mentioned above, a deterministic jitter ja 2 -Hi=0.2 Ul is also approximately recovered as 
it has been input for the present test PDF before. 

15 A second, refined embodiment of the invention deals with fitting an n-th order polynome 
to the BER. Doing this said small offset due to interpolating the quadratic function with 
its top touching a baseline B with a log-io(BER) level of 0 (probability of an event 

equals 1) according to the previous embodiment can be accounted for. The polynomial 
function to be fitted to the BER reads as: 

20 cp(x) = ag 4- a^x + . . . 4- a m x m . 

To determine the coefficients the method of least mean squares can be used. In order 
to keep the computational effort limited a restriction to a quadratic function seems to 
be sufficient: 

25 

(p(x) = + aix + a2X 2 

The linear equation system that needs to be solved takes the following form: 
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wherein the Pj are the measured data points used for the fit. The equation system can 
be solved according to Cramer's rule. 

5 According to this embodiment the trendline T is computed by applying a coordinate 
transform of the fitted quadratic function with the goal of translating the origin to the 
summit of the parabola that represents the second order polynomial fit function. 
Thereafter the same technique of deriving the random and deterministic jitter is 
deployed in the shifted coordinate system as it was the case in the simpler version of 
10 the previous embodiment. 

What is important for the measurement time is the fact that by applying the method of 
the present invention one is able to represent the random jitter characteristic by means 
of a quadratic approximation, while the method also allows an extrapolation in order to 
extract an amount of total jitter TJ at a given BER level. By just taking a few 

15 measurement data, e.g., at a BER of between 10-4...10- 6 , and then to extrapolate for 
the total jitter at a second baseline B-|2 having a BER-value of 10- 12 , measurement 
time is considerably reduced. An embodiment considering this approach is illustrated 
in Fig. 5. The polynomial fit of second order is applied to BER-measurements 
(reference numeral "M" in Fig. 5) taken at three timing points. No offset has been 

20 considered in this embodiment. The theoretical curve regarded as being "measured" 
in this document is once more displayed as a thin line for comparison. 
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